Control of heat transport in quantum spin systems 
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We study heat transport in quantum spin systems analytically and numerically. First, we demon- 
strate that heat current through a two-level quantum spin system can be modulated from zero to 
a finite value by tuning a magnetic field. Second, we show that a spin system, consisting of two 
dissimilar parts - one is gapped and the other is gapless, exhibits current rectification and negative 
differential thermal resistance. Possible experimental realizations by using molecular junctions or 
magnetic materials are discussed. 
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When an electron moves in solids, it carries charge, 
spin, as well as heat. The charge degree, the basis of the 
modern electronics, has been fully studied over the past 
decades and it is still being under intensive investigation 
at the molecular level [l|, |j|. On the other hand, the 
spin degree of freedom, which may be utilized to carry 
and process information as well, has also been explored 
intensively in past years, which has resulted in an emerg- 
ing field "spintronics" [3|ji4l- I'^ ^^^^ context, spin (mag- 
netization) current in insulating magnetic materials has 
attracted considerable attention [sj, |6[. However, heat 
current in spin systems has not been well studied as it 
should be. 

In fact, with the rapid miniaturization and increase 
of operational speed of microelectronic devices, a great 
amount of redundant heat are produced, which will in 
turn affect device performance. Thus heat dissipation 
and heat management is becoming more and more impor- 
tant [3,|8|. Moreover, it was found recently that heat due 
to phonons can be used to carry and process information 
[^]. Therefore, the study of heat conduction is not only 
of fundamental important but also helpful for the design 
and fabrication of heat dissipator and phononic devices. 
Indeed, many interesting conceptual devices such as ther- 
mal rectifier (diode) [10|, lul lid. | l3l| , thermal transistor 
lj,ll5|, and thermal logic gate [16| have been proposed. 
Experimentally, a nanoscale solid state thermal rectifier 
using deposited carbon nanotubes has been realized re- 
cently [l7| , and a heat transistor - control heat current of 
electrons- controlled by a voltage gate has also been re- 
ported [lg|. Furthermore, Segal and Nizan showed that 
thermal rectification can appear in a two-level system 
asymmetrically coupled to phonon baths [l9|, |20|, pro- 
viding the possibility to control heat at a microscopic 
level. 

In this Letter, we demonstrate that the heat part from 
spins can be also modulated and controlled. For example, 
we show that heat current in a two-level system can be 
modulated from zero to a finite value by tuning the mag- 




FIG. 1: (color online) A schematic representation of the model 
with size A'^ = 6. The spin model is connected to two phonon 
baths held at different temperatures, T^ and Tr. An inhomo- 
geneous field is appHed to introduce an asymmetric structure. 



netic field, h. Near h ~ Q, the two levels are almost de- 
generate and the system can jump easily between them; 
thus, heat current (proportional to h?) is small. As h in- 
creases from zero to a finite value, heat current increases 
accordingly. We further consider a spin-1/2 system con- 
sists of two different segments: one is gapped and the 
other is gapless. We show that in such a structure ther- 
mal rectifying efficiency can be very high, up to 10 times. 
This system also exhibits negative differential resistance, 
a feature which is necessary for building up a thermal 
transistor. 

Model For simplicity, we consider an inhomogeneous 
mesoscopic spin-1/2 chain whose Hamiltonian reads 
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where N is the number of spins, the operators a^ and 
a^ are the Pauli matrices for the nth spin, Q is the cou- 
pling constant between the nearest neighbor spins, and 
hn is the magnetic field strength (Zeeman splitting) at 
the nth site. Fig. [1] shows a schematic representation of 
this model. In the thermodynamic limit, the model is an 
analog of the Hcisenbcrg spin-1/2 chain, wherein the heat 
conduction is ballistic [2l|, [2^ . At the microscopic level, 
([1]) may be viewed as a sim plif ied version of a molecular 
model in low temperatures [20|. 

We use the quantum master equation (QME) to study 
heat conduction in this model. Two phononic baths of 
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FIG. 2: (color online) Modulation of Heat Current. Heat 
current as a function of the field. The bath temperatures are 
Tl = To + 0.05 and Tr = To - 0.05. Three cases are shown: 
To = 0.2 (solid line), To = 0.4 (dashed line), and Tb = 0.6 
(dotted line). 



different temperatures are connected to the system at the 
ends. By tracing out the baths within the Born-Markov 
approximation, we obtain the equation of motion for the 
reduced density matrix of the system iji — 1) |23l.l24|. 



dp 
'dt 



-i[H,p] +ClP + Crp, 



(2) 



where ClP {Crp) is a dissipative term due to the cou- 
pUng with the left (right) bath. ClP is given by ClP = 
{[XlP,cfi] -I-H.c.), and C^p can be given in a similar 
way. Here the operator Xl is defined through 

{■m\XL\n) = Xem,nnL{emji){ni\(jl\n) , (3) 

where £„,« = £m -£«, and nL{em,n) = (e^"-"/^^ - l) 
is the Bose distribution (fcs = 1) with Tl being the tem- 
perature of left bath. {|n)} and {en} are the eigenstates 
and eigenvalues of the system Hamiltonian 7f , respec- 
tively. A is the coupling strength with the bath. The bath 
spectral function we used is an Ohmic-type. The master 
equation ^ is solved by the fourth-order Runge-Kutta 
method. In numerical simulations, we take A = 0.01; the 
simulation time is chosen long enough such that the final 
density matrix p reaches a steady state. 

Modulation of heat current To demonstrate the modu- 
lation of heat current by an external field, we consider a 
simple two-level system (TLS), namely, iV = 1; thus Eq. 
([T]) becomes H = ha^- . From the QME, we obtain heat 
current. 
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Fig. [2] shows heat current as a function of the field. 
The temperatures of the left and right baths are Tl = 
To + 0.05 and Tr = Tq - 0.05, respectively. To is the 
mean temperature. In Fig. ^ we see that heat current 
first increases with the field and then decreases. In low 
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FIG. 3: (color online) Rectification of Heat Current, (a) 
Heat current vs temperature difference. The temperatures of 
the baths are Tl = To{l + A) and Tr = ro(l - A) where 
To is the mean temperature, ranging from To = 1 to To = 5. 
iV = 6 and Q = 0.6. /i„ = 1 ff 1 < n < N/2 and /i„ = 
otherwise, (b) Normalized heat current [J/J(A = 0.5)] vs 
A with different couplings. Here To = 1. The insets show 
the rectification efficiency, rj = |J+/J_|, vs the temperature 
difference. The lines are to guide the eye. 



fields (h <C To), ul^r ~ Tl^r/H^ then J oc h- . In high 
fields [h > To), UL^R '^ e-^''/^^.^, then J oc h^e-^''^'^°, 
implying that heat current decays to zero when h is large. 
Therefore, in such a model we can modulate the current 
from zero to a finite value by gradually switching on the 
external field. This result is similar to those observed 
in one dimensional spin- 1/2 systems recently [25l|26|, al- 
though in our case we consider just one spin. In Fig. [51 we 
can also see that the current can be modulated in a much 
wider region in a high temperature, such as Tq = 0.6. 

Rectification of Heat Current We now consider N spins 
in an inhomogeneous magnetic field, namely, /i„ = 1 (as 
the unit of energy) if 1 < n < N/2 and /i„ = oth- 
erwise (see Fig. [T]). The temperatures of the left and 
right baths are Tl = To(l -t- A) and Tr = To(l - A), 
respectively, where Tq is a mean temperature and A is 
the dimensionless temperature difference. The current 
operator may be defined through the equation of conti- 
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FIG. 4: (color online) (a) Energy gap of an isolated homo- 
geneous system (ft„ = h) vs the system size. The gap is the 
energy difference between the ground state and the first ex- 
cited state. Q = 0.6. (b) The probability, Po, to find the left 
part of the system in the ground state as a function of A. 
To = 1, iV = 6, Q = 0.6, and /i„ = 1 if 1 < n < N/2 and zero 
otherwise. 



niiity. In Fig. [31 we show the heat current as a function 
of the temperature difference with the mean temperature 
ranging from Tq = 1 to Tq = 5. When the temperature 
is low (To = 1), we observe that for A > the heat cur- 
rent increases with A, while in the region A < the heat 
current remains very small. Thus, our model exhibits 
thermal rectification; namely, heat flows favorably in one 
direction. However, when the temperature becomes high, 
such as To — 5, the magnitude of heat current changes 
little as the bath temperatures are exchanged. In this 
case the model cannot act as a good rectifier. Neverthe- 
less, in a wide range of temperature (Tq < 3), this model 
shows thermal rectifying effect; the mechanism will be 
illustrated later. In Fig. El^b), we show the heat cur- 
rent for a model with different couplings. We can see 
that the rectifying effect may sustain to large coupling 
constants. To quantify the rectification efficiency, we in- 
troduce the ratio, rj = \J+/J^\, where J+ is the current 
when Ti > Tr and J_ is the current when the temper- 
atures are swapped, i.e., Tl < Tr. In a weak coupling 
case, e.g., Q = 0.2, the efficiency may be more than 10 
[see insets of Fig.[21[a) and (b)]; however, as the couphng 
Q becomes stronger, the efficiency decreases. 

To understand the rectifying effect qualitatively, we 
calculate the energy gap Eg for a homogeneous isolated 
system (/i„ = h) and show it in Fig. HJa). The energy 
gap is obtained by calculating the energy difference be- 
tween the ground state and the first excited state. In a 
finite field ft, = 1, we observe an energy gap, Eg r^ 2.8, 
which is slightly dependent on the size. In the absence 
of the field {h — 0), a small gap which vanishes in the 
thermodynamic limit is seen. Note that if the size iV is 
odd, there is no gap (degeneracy of the ground state); 
this is due to the particle-hole symmetry if we map the 
spin model to a free Fermion one by the Jordan- Wigner 
transformation. 

At low temperatures, only a few low-lying states of 
each part are relevant. Therefore, if the left part (with 
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FIG. 5: (color online) Rectification of Heat Current. 

Heat current of a system with difference sizes. The bath tem- 
peratures are Tl = To(l + A) and Tr = To{l - A). To = 1 
and Q = 0.6. The field is /i„ = 1 if n < N/2 and /i„ = 
otherwise. The inset shows the rectification efficiency as a 
function of the temperature difference. 



a gap) is in contact with a cold bath and the right with 
a hot bath, i.e., T^ = Tc and Tr — Th [Tc < T/i), then 
the left part mainly remains in the ground state. This 
is reflected in Fig. IJ^b), which shows the probability, Pq, 
to find the left part in the ground state for different A. 
In this case, the transition rate of the left part between 
different levels is low, and thus the heat current is small 
although the right part (no gap or a small gap) may jump 
easily between different levels (see also Fig. [2] for refer- 
ence) . Reversely, if the left part of the system is in con- 
tact with a hot bath and the right with a cold one, the 
transition rate of the left part between different levels 
is large, and then heat current becomes large. This can 
also explain the low rectifying efficiency when the mean 
temperature Tg is increased. In this case, the transition 
probability of the left part between different levels may 
also be large; then, the magnitude of heat current changes 
little when the bath temperatures are exchanged, imply- 
ing a low efficiency. In fact, we may observe thermal 
rectification provided that Tq ^ Eg. 

In Fig. [5l we show the heat current for a system with 
different sizes. Here the coupling Q — 0.6, and the field 
is /i„ = 1 if n < N/2 and hn — otherwise. In the small 
size case (N = 4), we just observe thermal rectification 
with a low efficiency. The reason may be that the effec- 
tive coupling between the left and the right parts can be 
strong for a small size system. As a result, the left part 
may be excited by the right part even though it is con- 
nected to a cold bath. However, in a larger size system, 
i.e., A^ = 6 or iV = 8, we may observe both thermal rec- 
tification and negative differential resistance. Note also 
that the efficiency changes very little when A^ = 6 or 
N ^ 8, implying that the model may act as a rectifier in 
an even larger size case. 

Negative Differential Thermal Resistance (NDTR) In 
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FIG. 6: (color online) Negative differential thermal re- 
sistance (NDTR). Heat current vs temperature of left heat 
bath, Tl, with fixed Tr — f.5. The other parameters are 
Q = 0.6, Af = 6, /i„ = 1 if 1 < n < N/2 and /i„ = other- 
wise. The NTDR is clearly seen for Tl < 1.0. 



fact, in Fig. [31 we can also observe NDTR in the re- 
gion of A < 0, i.e., the decrease of heat current with the 
increase of temperature difference. A clearer representa- 
tion is shown in Fig. [6l where we fix the temperature of 
the right bath, Tn = 1-5. We see that when the tem- 
perature of the left bath Tl is increased from 0.4 to 1.0, 
i.e., decreasing the temperature difference, thermal cur- 
rent increases. The reason is that: If Tl is low, the left 
part is rarely excited, implying a small current; other- 
wise, current may become large. NDTR is an important 
physical property that may be used to build spin-based 
thermal transistors. 

Summary We have studied the possibilities to control 
heat current in mesoscopic spin models. We have showed 
that heat current could be modulated from zero to a fi- 
nite value in a two-level system by tuning the magnetic 
field. We have also studied thermal rectification and 
negative differential thermal resistance in an asymmet- 
ric model. The model consists of two parts: the left part 
is gapped, and the right part is gapless. Such a structure 
is of great importance for the model to exhibit rectifica- 
tion and NDTR. In certain cases, we have found that the 
rectification efficiency, |J+/J_|, can be larger than 10. 
Finally, we would like to discuss the possible realizations 
of the model in experiment. The first is to make use of 
the asymmetric structure in molecular bridges that can 
be easily introduced. For example, we may use a molecule 
consisting of two (weakly) coupled nonidentical spatially 
separated segments; each is taken to be an anharmonic 
system, e.g., anharmonic vibrations or molecular libra- 
tions, where at low temperatures only the lowest (two) 
quantum states are relevant (see Fig.[T|). However, in this 
case, there are gaps in both segments, so the rectifying 
effect may be not so high. The second possible way is to 
use magnetic materials or molecular magnets 27[. The 
model may be made up of two coupled magnetic mate- 
rials: one is gapped and the other is gapless. For the 



gapped material, one could use spin-Peierls systems or 
introduce a magnetic field to open a gap {h/Q > 2) [2q . 
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